Depletion forces in non-equilibrium 
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The concept of effective depletion forces between two fixed big colloidal particles in a bath of 
small particles is generalized to a non-equilibrium situation where the bath of small Brownian 
particles is flowing around the big particles with a prescribed velocity. In striking contrast to the 
equilibrium case, the non-equilibrium forces violate Newton's third law, are non-conservative and 
strongly anisotropic, featuring both strong attractive and repulsive domains. 
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Colloidal mixtures are excellent model systems to 
study fundamental questions of equilibrium and non- 
equilibrium phase transitions with important physical 
applications, such as protein crystallization 1] as well 
as gel and glass formation Q • One of the reasons for the 
pivotal role of colloids as model systems lies in the fact 
that their interactions can be tailored Q . In the theoret- 
ical description of colloidal mixtures, the concept of the 
effective interaction has proved to be of paramount im- 
portance in advancing our understanding of how to steer 
the properties of these complex systems :4j. It stems 
from integrating out the coordinates of the small parti- 
cles, reducing thereby the number of statistical degrees 
of freedom considerably 0, |(| . This approach has been 
successfully applied to colloid / polymer mixtures 0,13, @> 
to binary colloidal mixtures \WL Hll ] and to charged sus- 
pensions [l^. The most widely known effect is the de- 
pletion attraction in strongly non-additive systems, such 
as the Asakura-Oosawa model Q pertaining to colloids 
and non-adsorbing ideal polymer coils. This attraction is 
entropic, i.e., it scales with the thermal energy fc^T and 
results from a depletion zone of polymer between a pair 
of nearly touching colloidal particles, which causes an 
imbalance of the osmotic pressure of the polymers acting 
on the colloidal surface. While the range of the depletion 
attraction is typically associated with the radius of gyra- 
tion of the polymers, its attractive depth is fixed by the 
polymer concentration, so that both the range and the 
depth are in principle tunable. 

The crucial limitation of the concept of effective inter- 
actions is that it works, strictly speaking, only in thermo- 
dynamic equilibrium and can thereafter be applied for the 
calculations of static quantities, such as correlation func- 
tions and free energies (resulting into phase diagrams). 
In this Letter, we investigate whether and how the ef- 
fective interaction can be applied also in non- equilibrium 
steady state situations. In order to be as transparent as 
possible, we study a minimal model of non-equilibrium 
for two fixed big colloidal particles in a quiescent solvent 
which are exposed to a flowing bath of small Brownian 
particles. The latter ones arc ideal, in full analogy to 
to the equilibrium case of the Asakura-Oosawa model, 
and hydrodynamic interactions mediated by the solvent 



are neglected. Despite its simplicity, the model is real- 
ized in colloidal mixtures with weakly interacting dcple- 
tants (which justifies the assumed ideality) which have 
a small physical volume fraction with well-separated ex- 
cluded volumes of the big and small particles (which jus- 
tifies neglecting hydrodynamic interactions). Examples 
include strongly asymmetric highly c harg ed suspensions 
with a small physical volume fraction 13J and sterically- 
stabilized colloidal particles in a bath of thin rods or 
polymer coils 0, . Our situation is different to the 
classical hydrodynamic problem of two fixed spheres in a 
Stokes flow of a viscous incompressible fluid, |17l Ull fl9| 
in which the force vanishes for zero flow. More impor- 
tantly, the physical origin of the depletion forces is the 
high osmotic compressibility of the bath particles com- 
plementary to the incompressibility of the viscous fluid 
mediating hydrodynamic interactions. 

Based on a standard Smoluchowski approach and 
on Brownian dynamics computer simulations in non- 
equilibrium, we show that the effective force can be gen- 
eralized in a straightforward way to a stationary non- 
equilibrium state. The results for the effective force, 
however, differ from the corresponding equilibrium case 
both qualitatively and quantitatively. First of all, the 
non-equilibrium force field (as a function of the mu- 
tual distance between the particles) is non-conservative 
and violates Newton's third law, implying that the con- 
cept of an effective interaction potential does not hold 
in non-equilibrium. Second, when the small particles are 
brought to flow, the equilibrium depletion attraction gets 
strongly anisotropic favoring an alignment of the colloidal 
pair along the flow direction. Conversely, there is mutual 
effective repulsion between the colloidal pair perpendicu- 
lar to the drift direction of the bath due to a compression 
of the small particles between the colloids. Furthermore, 
the range of the depletion interaction in the drifting case 
is larger than that in equilibrium. 

Equilibrium depletion forces can be measured in real 
colloidal samples with high accuracy 0, 0, 0, |2(3, . 
Our predictions in non-equilibrium can be tested in ex- 
periments of binary colloidal mixtures as well. A drifting 
bath is generated either by gravity [U H3| or by an elec- 
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FIG. 1: Steady-state contour density fields of non- 
interacting Brownian particles, driven around two sta- 
tionary colloidal spheres, on the x — z plane where 
the centers of the spheres are located. The distri- 
butions are obtained from the superposition approxi- 
mation, Eq. {T|l, Bright regions show high densities 
while low density regions are plotted dark. One col- 
loid is placed at the origin, Ri = 0, and the posi- 
tion R.2 of the other one is: (a) R 2 / f = (2,0,0); 
(b) R 2 /a = (1.5,0,1); (c) R 2 /a = (1,0,2); (b) 
R2/C — (0,0,4). In this figure, the bath particles drift- 
ing from the left to the right with a speed v* = 5. 



fixed by optical tweezers pH |24| . The measurement of 
the relative displacement with respect to the tweezer's 
position of the big particles provides direct access to the 
effective non-equilibrium forces. 

In our theoretical model we consider two fixed big par- 
ticles in a bath of small Brownian particles that drift with 
velocity v = vz. The bath has a bulk number density 
Po . While there is no interaction between the bath parti- 
cles, the interaction between the big and small particles 
is described with an exponential, spherically symmetric 
potential V(r) given by 

V(r) = V(r)=V e X p[-(r/ ( r) n } ) 

where r is the vector separating the bath particle with 
the center of the colloid, r is the magnitude of r and 
Vq = lOksT. In this work we focus on an exponent 
n = 6, on one hand to deal with a well-defined 'hard' 
interaction diameter of the colloidal sphere, a, on the 
other hand to circumvent numerical inconveniences due 
to huge forces at r = a for high exponents n. 

In what follows, we first discuss a theoretical Smolu- 
chowski approach for a single driven colloid. Adopting a 



FIG. 2: The depletion force F<j e pi,2 (arrows) acting on 
colloid 2, when colloid 1 is placed at the origin (dashed 
circle), as the Brownian bath particles drift from the left 
to the right with speed v* . Forces are plotted in the 
X — Z plane where the centers of the colloids are lo- 
cated. Results are shown for drift velocities (a) v* = 
(equilibrium), (b) v* — 1, and (c) v* — 5. The results in 
(a)-(c) are calculated theoretically from Eqs. JIJ and J2J). 
Plot (d) shows the force for v* = 5 resulting from Brow- 
nian Dynamics computer simulations. Different scales 
for the length of the vectors are used for different val- 
ues of v*. The magnitude of the force vector located 
at (X,Z) = (0,2) is l3Fo = 1.22 in (a), 2.38 in (b), 
8.04 in (c), and 9.16 in (d). Furthermore, po = l/cr 3 . 



dynamical superposition approximation (DSA), we then 
address the case of two colloids and afterwards we com- 
pare the theoretical predictions arising from the DSA 
with Brownian dynamics (BD) computer simulations, 
finding good agreement between the two approaches. In 
detail, for the case of a single big colloidal particle, the 
Smoluchowski equation [19| for the steady state density 
field p(r) of the bath particles reads as V ■ j(r) = with 
the current density field j (r) composed of 

j(r) = - Vp(r) + I>(r)W(r) + vp(r), 

where /3 = 1 /ksT is the inverse thermal energy and T the 
mobility of the bath particles in the solvent. As the prob- 
lem possesses azimuthal symmetry, these equations can 
be readily numerically solved in cylindrical coordinates. 
Let us define the resulting density profile of the bath 
particles around one big colloid as p^'(r — R; v) with R 
denoting the position of the colloid. Spatial homogeneity 
dictates that p^ depends only on the difference r — R. 
The force Fd acting on the big particle due to the drifting 
bath is then obtained as Fd = — j d 3 rp^ (r)W(\r\). 

In the presence of two colloids, the correspond- 
ing quantity of interest is the bath density profile 



p^(r, Ri, R2; v) that depends on the locations R; of 
the colloids, (i — 1,2), and parametrically on the ve- 
locity v. In this case, the numerical solution of the 
Smoluchowski equation is more difficult such that a di- 
rect particle-resolved computer simulations is more ap- 
propriate. However, let us motivate a simple superposi- 
tion approximation. As the bath particles are noninter- 
acting, in the case of thermodynamic equilibrium, v = 0, 
p^(r, Ri, R 2 ; 0) can be exactly factorized as the prod- 
uct of the density fields p^(r — R 4 ; 0) divided by p . In 
what follows, we employ the same factorization ansatz 
also for the non-equilibrium situation, v 7^ 0. Thereby, 
we introduce the dynamical analog of the superposition 
approximation that was first introduced by Attard |2^] 
for equilibrium problems dealing with interacting parti- 
cles, in order to calculate static depletion forces. Hence, 
in this dynamical superposition approximation (DSA) , we 
write: 

pW (r, R l5 R 2 ; v) « p {1) (r - R l5 v)p« (r - R 2 ; v)/p .{\) 

Results from the DSA are shown in Fig. Q where we 
plot the density distribution in the x-z plane, on which 
the centers of both colloids are located. The profiles are 
plotted for four different relative positions of the two col- 
loids and for a dimensionless drift velocity v* = 5, de- 
fined as v* = (3o~v/T. The density profiles are strongly 
anisotropic and show long-ranged low-density regions in 
the wake region behind the spheres. The extension of 
the concept of the effective interaction to a steady-state 
non-equilibrium situation can now be performed via the 
definition of the effective force which is a statistical av- 
erage over the forces of the small bath particles exerted 
onto the big ones p|- Generalizing to a steady-state av- 
erage in non-equilibrium, the effective force reads as 26]: 

F i = ~[ d 3 rp( 2 )(r,R 1 ,R 2 ;v)V Rl V(|R l -r|). (2) 

The force has been calculated in theory by employ- 
ing the DSA, Eq. ([[J, for various different configura- 
tions of the colloids and drift velocities, including the 
equilibrium case, v — 0, in which the results are ex- 
act. In this way, a comparison between the static and 
dynamic depletion force F<iepl,i acting on the z-th parti- 
cle can be obtained. The latter is defined as the force 
Fdepi,i = Pi — Fd, relative to the constant the drift force 
Fd acting on a single big particle (fiF^a = 2.96 for v* = 1 
and /3-FdC = 11.76 for v* — 5). We accompanied the 
DSA with standard Brownian Dynamic simulations |27| , 
in which the Langevin equations for each bath particle 
are solved iteratively. Here, the bath particles are driven 
by an external force v/T for a prescribed flow velocity v. 
The box lengths in all simulations were L x = L y = 10a 
and L z = 20a and N = 2000 particles were simulated, 
applying periodic boundary conditions in all directions. 

Results for the averaged force field acting on a colloid 
placed at the position R 2 away from another colloid at 
the origin (Ri = 0), are depicted in Fig. El A comparison 
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FIG. 3: Depletion forces as in Fig. but now we 
plot the projection of the the depletion force F*(X) = 
/3crFdo P i,2 • X for fixed Z = in the X direction (a) 
and the projection of the depletion force F*(Z) = 
/3crFdcpi,2 • z for fixed X = in Z-direction. The solid 
line is the exact result for v* = and the symbols 
are results from the DSA for v* = 1 (squares) and 
v* = 5 (diamonds), as well as simulation results for 
v* = 1 (triangles up) and v* = 5 (triangles down). 



between Figs. I^c) and[2{d), (DSA vs. BD-simulation) , 
shows that the former captures the features of the non- 
equilibrium depletion force quantitatively, even for drift- 
ing speeds as large as v* — 5. Whereas for the case 
of equilibrium, Fig. |2Ja), we obtain spherically sym- 
metric attractive depletion forces, the dynamical deple- 
tion force shows new qualitative features: First, as is 
clear from inspection of Fig. [5J the dynamical deple- 
tion forces violate Newton's third law. Second, the vi- 
olation of the "action = reaction" -principle implies a 
non-conservative force. Indeed, if a putative depletion 
potential Vd cp i existed in the dynamical case at hand, 
translational invariance of space (which is not broken by 
the steady flow) would imply Vd ep i = Vd op i(Ri — R2)- As 
FdcpM = -V Ri Vde P i(Ri - R2), i = 1, 2, this would have 
the consequence Fd op i,i = — Fd cp i.2- Third, the deple- 
tion force shows strong anisotropy, becoming attractive 
as the colloid starts positioning itself 'behind' the one at 
the origin and its magnitude (in the attractive case) and 
its range are vastly larger for v ^ than in equilibrium. 

In order to better quantify the features of the deple- 
tion force in non-equilibrium, we show in Fig. 3 plots 
of the projected depletion force at fixed X = (parti- 
cles fully aligned along the flow direction) and at fixed 
Z = (particles fully aligned perpendicular to the flow 
direction). The depletion force exhibits a repulsive bar- 
rier in the A-direction due to compression of the small 
particles between the colloids, in qualitative difference 
to the purely attractive equilibrium case v = 0. In the 
direction, on the other hand, there is much stronger 
attraction for Z > due to a long low density region in 
the wake of the colloids but a marked repulsive barrier for 
Z < 0. Furthermore, the range of the depletion force is 
significantly larger than in equilibrium. Notice that we 
show forces also for the case of completely overlapping 
colloids [13, -2.0 < Xj a < 2.0 and -2.0 < Zja < 2.0. 
For Ri = R 2 , the depletion force on both colloids is ob- 
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viously the same, Fdcpi,i = Fd cp i.2 but, in contrast to 
the equilibrium case, it does not vanish. We can finally 
deduce from Fig. that the computer simulation results 
are in good quantitative agreement with the theory up to 
drifting speeds of v* — 1, but for large drifting speeds as 
high as v* = 5 deviations near the maxima and minima 
are getting more pronounced. 

In conclusion, we have generalized the depletion con- 
cept to a simple non-equilibrium situation of two big par- 
ticles in a Brownian flow of small ideal particles. The ef- 
fective forces averaged in the steady state are strikingly 
different from those in equilibrium: they are strongly 
anisotropic with both attractive and repulsive domains 
and violate Newton's third law, in strong resemblance to 
the so-called "social forces" that have been applied ex- 
tensively in the modelling of pedestrian and automotive 
traffic flows p8|. 

Although the present study has focused on noninter- 
acting bath particles, the extension to interacting ones 



can be carried out within the framework of the recently 
proposed dynamical density functional theory (DDFT) 
[29| . We expect that interaction effects will weaken the 
strength of the dynamical depletion forces, as they will 
tend to smoothen out regions of spatial inhomogencity 
of the bath density that give rise to these forces in the 
first place. Finally, we remark that the strong attrac- 
tive forces tending to align the colloids in the direction 
of the incoming flow can be regarded as precursors of 
the laning transition that has been analyzed recently for 
concentrated colloidal mixtures |30| . 
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